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RESPONSE ENVELOPE - A GLOBAL DESCRIPTION OF THREE-AXIS 


LARGE-ANGLE SPACECRAFT ATTITUDE CONTROL SYSTEMS 
By George Meyer 
Ames Research Center 


SUMMARY 


Arbitrary three-axis, large-angle attitude control systems of highly 
maneuverable spacecraft capable of tracking time-varying target attitudes are 
considered. Such systems are inherently multivariable, multidimensional, non- 
linear, and nonautonomous with infinite sets of forcing functions and initial 
conditions. The concept of a response envelope is introduced. It provides an 
analytical procedure for judging whether a proposed attitude control system of 
this type is fast enough for the mission of the spacecraft. It is defined as 
a scalar function of time which at each instant is the maximum possible three- 
dimensional attitude error between spacecraft and target for any admissible 
initial condition, time variation of target attitude, and disturbance. A sim- 
ple procedure is presented for computing the response envelope approximately. 
The approximation is conservative and may be used as a basis for accepting a 
proposed system. The procedure is based on the Liapunov theory and a proto- 
type of general attitude control systems. The proposed procedure is 
illustrated by a numerical example based on an Orbiting Astronomical 
Observatory. 


INTRODUCTION 


Systems for controlling spacecraft attitude vary greatly in their inter- 
nal structure. Thus, for example, torque may be generated by reaction wheels, 
control moment gyros, reaction jets, or by some interaction with the environ- 
ment (i.e., gravity gradient, solar pressure, magnetic field, etc.). Simi- 
larly, spacecraft attitude may be measured with star trackers, sun sensors, 
inertial gyros, or by using the earth’s magnetic field. Spacecraft angular 
velocity may be measured directly, or it may be computed from attitude data. 
Finally, the feedback linking the outputs of attitude and angular velocity 
sensors with the inputs to the torquers may be designed using Euler angles, 
some property of sensors (see refs. 1 and 2), or Euler’s theorem on three- 
dimensional rotations (see refs. 3 and 4). The design may be optimal in some 
sense, or it may simply be intuitively appealing. However, despite this 
diversity, every attitude control system is built for the single purpose of 
allowing the spacecraft to fulfill its mission. Consequently, there is always 
a stage in the design of such systems at which one must judge the quality of 
the proposed system relative to the mission of the spacecraft. Since all pos- 
sible situations must be considered, quality must be judged on the basis of 
the overall (global) properties of the system. This poses no difficulty if 
the mission of the spacecraft presents its control system with only a small 



number of situations to control because then global behavior may be obtained 
by directly testing the system or its analytical model. If, however, the num- 
ber of situations to control is so large that direct enumeration is impracti- 
cal, or, indeed, impossible, tests consisting of a small sample may not be 
decisive: there is no certainty that every case which results in mission 

failure is included in the test sample. 

The present note is concerned with the responsiveness of three-axis, 
large-angle attitude control systems of highly maneuverable spacecraft capable 
of tracking time-varying target attitudes. A system of this type is inher- 
ently multivariable, multidimensional, nonlinear, and nonautonomous with infi- 
nite sets of forcing functions and initial conditions. Since, on the one hand, 
there are infinitely many cases involved, and, on the other hand, there are no 
typical cases (because of nonlinearity) from which to extrapolate to any other 
case, the responsiveness of such systems cannot be determined by direct enu- 
meration of cases. Another method must be found. The purpose of the present 
report is to present such a method. 

The method involves two ideas, namely, the error angle first introduced 
in reference 4, and the response envelope introduced in the present note. The 
error angle is a distance between three-dimensional rotations. The response 
envelope is a function of time which at each instant is the maximum of all pos- 
sible values of the error angle at that instant. The response envelope indi- 
cates the responsiveness of the system in the following sense. Under no 
possible circumstances is the attitude error between spacecraft and target, at 
any instant, greater than the value of the response envelope at that instant. 

The discussion proceeds from the general to the particular in three steps. 
First, the response envelope is defined precisely. Second, a procedure for 
computing a useful approximation of the response envelope is described. Third, 
a numerical example is presented as an illustration of the method. One of the 
control laws introduced in reference 4 is assumed for an Orbiting Astronomical 
Observatory, and the following questions are considered: (1) What is the 

responsiveness of the system to step changes in the target attitude, and how 
well does the system follow time-varying target attitudes? (2) How signifi- 
cant is gyroscopic coupling? How sensitive is the system to (3) external 
torque disturbances, (4) variations in system parameters, (5) changes in the 
form of the control law, and (6) time delays in the controller? Answers are 
obtained by the proposed procedure. 


SYMBOLS 


A output matrix; actual attitude of spacecraft relative to inertial 

space 

A^ s input matrix; desired attitude of spacecraft relative to inertial 
space 

a orthonormal triplet of vectors fixed to the spacecraft 
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Ui (t) 
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unit eigenvector of R; error axis 

orthonormal triplet of vectors fixed to the target 

time derivative along a system trajectory 

set of system state equations (i.e., x = f [x, Ui (t) ,u 2 (t) ] ) 
right-hand side of system state equations 
scalar functions appearing in equation (16) 
angular momentum storage capacity of controlling device 
unit matrix 

spacecraft coordinates of moment of inertia of main body 

maximum eigenvalue of J a 

minimum eigenvalue of J a 

torque capacity of the controlling device 

set of all motions of the system 

perturbation function 

error matrix defined by equation (2) 

set of all real n-tuples 

matrix function defined by equation (A2) in appendix A 
orthonormal triplet of vectors fixed in inertial space 
time 

trace of matrix ( ) 

set of admissible time variations of u. 

—i 

value of u^ at t 

solution of Hamilton-Jacoby equation 
solution of Liapunov inequality 
explicit part of V + (t 5 x) 
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set of admissible time histories of spacecraft angular velocity 
body coordinates of inertial angular velocity of spacecraft 
spherical bound on w a 

target coordinates of inertial angular velocity of target 

spherical bound on w d 

state space of the system 

element of state space 

control law given by equation (16) 

region of operation of the system 

real constant 

set of admissible time histories of spacecraft acceleration 
angular acceleration of the spacecraft 
spherical bound on a a 
error angle 

point on a time history of error angle 
value of the response envelope at t 

value of lower estimate (of the response envelope) at t 
value of upper estimate (of the response envelope) at t 
time history (thus, x - {[x(t),t]: t >_ 0}) 

time derivative of ( ) 
matrix transpose of ( ) 


DEFINITION OF THE RESPONSE ENVELOPE 


Consider in general terms a complete attitude control system. It con- 
sists of essentially three distinct parts: (1) an input generator, (2) a 

spacecraft, and (3) a disturbance generator. The input generator represents 
all admissible time variations of attitude to be followed by the spacecraft. 
The spacecraft consists of the spacecraft inertia, torquers, sensors, and a 
controller which sends commands to the torquers based on information supplied 
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by the sensors. The disturbance generator represents all admissible time vari- 
ations of undesirable but unavoidable disturbances acting on the spacecraft 
such as perturbations in system parameters and external torque disturbances. 
Thus, an attitude control system as defined in the present report is a com- 
plete spacecraft and its environment. The mathematical model to be assumed 
for such a system consists of the following items: (1) an n-dimensional 

state space X; (2) a region of operation 6 in X to which the motion of 
the system is restricted by physical limitations of system components; (3) two 
sets U_i and U 2 of forcing functions, where for i = 1, 2, lb is a set of 
piecewise continuous vector functions of time uu with values u^(t) in 
Ui(t); and [4) a set E of state equations describing the dynamics of input 
generator, spacecraft, and output generator: 


E = {(Xi = fi[x,u 1 (t),u 2 (t)], i): i = 1, . . n} 

A motion of the system is a solution of E for some initial state x in 0, 
some input attitude generated by a forcing function ui in Uj , and some dis- 
turbance generated by a forcing function u 2 in U 2 . Thus, the sets E and 
0 x U ! x U_ 2 (Cartesian product) define the set M of all possible motions of 
the system. The central topic of this report is the description of M. Those 
aspects of M will be considered which reflect how well the spacecraft fol- 
lows any admissible time variation of the desired attitude in the presence 
of any possible disturbance. 

To decide at any instant of time how near the actual attitude of the 
spacecraft is to the desired attitude, it is necessary to have a notion of a 
distance between three-dimensional rotations. This may be introduced as fol- 
lows^ Consider three right-handed orthonormal triplets of vectors, say, s, a, 
and d. The attitude of a relative to s is given by the 3x3 orthogonal 
matrix A as whose elements are the direction cosines of a relative to s . 
Similarly, the attitude of d relative to s is given by A^ s . The attitude 
of a relative to d is given by the matrix A as A^j s , where t denotes 
matrix transpose. It follows (see ref. 5) from Euler's theorem on rotations 
that ma y he interpreted as a rotation about a single axis which is 

the eigenvector of A as A^ s . The angle of this rotation (according to para- 
graph 3 in appendix A) is the following function of A as and A^ s : 

'KAds^as) = arc cos ( y [tr (A as Ad s j - l] \ (1) 

[0,u] ^ J 

It is shown in appendix A that this function is a metric on the space of three- 
dimensional rotations. Indeed, (i) <|> (A<is > A as ) >_ 0; (ii) (j) (Ads > A as ) = 0 if 
and only if A as = Ad s ; (iii) <HAds,A as ) = <f> (Aas , Ads) ; Civ) if B is any 
three-dimensional rotation, then <f>(Ads> B ) + ‘I’CBjAag) >_ <KAd s , A as ) . Conse- 
quently, <p (Ads > A as ) may be interpreted as representing the distance between a 
and d. In fact, the following interpretations will be made in the remainder 
of this report. In inertial space s is fixed. In the spacecraft a is 
fixed; A as is the actual attitude of the spacecraft, and is the output of the 
attitude control system; Ad s defines the desired attitude, and is the input 
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to the system. The attitude error is the 3x3 orthogonal matrix R defined 
by 

R = A as Ads (2) 

It may be noted that R locates a relative to d. The magnitude of atti- 
tude error is * the error angle 4> = <J> (A^ s , A as ) . The direction of attitude 
error is the axis c (unit eigenvector) of R computed from R using 
equation (A5) in appendix A. 

It can be shown (see appendix A) that <J> has the following intuitively 
appealing properties. (1) It is the shortest unrestricted angular distance 
between the actual and the desired attitude. (2) It^is no smaller than the 
angle between the ith (i = 1, 2, 3) vector of the a-triplet and the ith 
vector of the 3-triplet. (3) In the one-dimensional case (i.e., shaft- 
positioning servos) in which rotation about only a single axis is allowed, the 

usual definition of error is <j>e = <f>d ~ where <f>d and <j>a are the input 

and output angles, respectively. In that case <J> = |<j> e | if |<j) e | <. tt, and 
4> = 2 tt - |cf) e | if tt <_ \$ e \ £_ 2tt. (4) When <J> is small, the attitude error 

may be represented by the vector <j>c whose components are the Euler angles of 

R, and whose magnitude <j> is the square root of the sum of the squares of 

these angles. For these reasons it appears that the error angle <j> is both a 

mathematically convenient and intuitively appealing scalar representation of 
three-dimensional attitude error of a spacecraft at each instant of time. 

The response envelope can now be defined in terms of time variation of 
error angle. With each admissible motion of the system it is possible to asso- 
ciate a time history of the error angle. For a given admissible initial 
state, and forcing functions u_i and u 2 , is a curve in the first quadrant 
of the t-cf) plane. The (^-coordinate of 4>_ will be denoted by <J> (t , x,^ ,u 2 ) . 
Let $ be the set of all such curves generated by M. The response envelope, 
to be denoted by $**, is defined to be the curve in the first quadrant of the 

t-<j> plane such that every point <J>**(t) of 4)** is the maximum of all values 

of 4^ in at time t. For this definition to make sense formally, it will 
be assumed that for all t >_ 0, the set 0 x U]_(t) x U 2 (t) is compact and that 
f [x,u 2 (t) ,u 2 (t) ] is uniformly bounded on this set. This assumption is suffi- 
cient to guarantee the existence of the response envelope, and represents an 
insignificant physical restriction. The value of the response envelope may be 
computed for any t >_ 0 by the following formula. 


4>**(t) = max [ max $ (t ,x,uj ,u 2 ) ] (3) 

x e 0 (u i , u 2 ) c Uj x U 2 

The following diagram summarizes the discussion of the present section. 
The input and disturbance generators generate admissible time variations of 
desired attitude Aj s and disturbance d_, respectively. Both act on the 
spacecraft to generate the rotation A as . For the complete system the follow- 
ing relations hold: the state space X C Xj x X 2 x X 3 , the region of 
operation 0 C e i x e 2 x 0 3 > an ^ the state equations E C Ej U E 2 U E 3 . By 
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INPUT SPACE- DISTURBANCE 

GENERATOR CRAFT GENERATOR 



Figure 1.- Attitude control system and its 
response envelope. 


means of equation (1), the distance between 
the desired and actual attitudes is computed 
to obtain the curve The boundary of all 

such curves is the response envelope jr**. 
Thus, the response envelope is a global prop- 
erty of an attitude control system that 
describes the responsiveness of the space- 
craft. Regardless of where in its region of 
operation the system is initially and regard- 
less of how the input attitude and distur- 
bance vary in time within the prescribed 
limits, the attitude error will not at any 
time be greater than the value of the 
response envelope at that time. 


APPROXIMATE COMPUTATION OF RESPONSE ENVELOPE 


Consider, now, the problems involved in the computation of the response 
envelope whose points are defined by equation (3) . The usual functional maxi- 
mization techniques yield in general only local maxima. In the absence of the 

a priori knowledge_of the number of these maxima, such methods will yield only 
lower estimates of the response envelope <£** , where <f) (t) <_ <j)**(t) for 

all t >_ 0. This is useful because an attitude control system may be judged 
unacceptable on the basis of such an estimate: there are conditions for which 

attitude error is at least as large as <f> (t) . In the present section a 
method is presented for computing an upper estimate <£_ such that 

<J>**(t) <_ cf>+(t) for all t >_ 0. This is useful because” an attitude control sys- 
tem may be judged acceptable on the basis of such an estimate: under no cir- 

cumstances can the attitude error be larger at time t than <)> + (t). The 
discussion proceeds in three steps. First, a general approximation procedure 
is set up using a state space interpretation of equation (3). Next, a proto- 
type of arbitrary attitude control systems is formulated. Finally, the compu- 
tation of an upper estimate of the prototype response envelope is stated 
explicitly . 


State Space Interpretation of the Response Envelope and a 
General Procedure for Computing Its Upper Estimate 

Let the multidimensional state space X and region of operation 6 be 
represented schematically as shown in the following figure, and let one of the 
state coordinates be the error angle <j> (transforming the state space if nec- 
essary) . The motion of the system starting in a state x and forced by uq 
and u 2 may be represented by a trajectory in 0, whose ^-coordinate at each 
instant is the value <J> (t ,x,uj ,u 2 ) of the error angle at that instant. The 
same initial state but a different pair of forcing functions results in a dif- 
ferent trajectory. Consider the bundle of such trajectories all emanating 
from the same initial state and generated by all admissible Uj and u 2 . This 
bundle defines a moving set of states which is reachable from the given state 
at a given time. The crosshat ched regions in figure 2 show schematically such 
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Figure 2.- 


Motion of* the set reachable 
from x. 


a set at various times. The maximum projec- 
tion 4>*Ct,x) of this set on the <j)-axis 
gives at each t the value of the bracketed 
term in equation (3). For example, in fig- 
ure 2, c|>*(t 2 ,x) is the maximum projection of 
the moving set at t = t 2 . Different ini- 
tial conditions generate different moving 
sets. Consider at each instant the union of 
all such moving sets generated by all ini- 
tial states in 0. This union defines a mov- 
ing cloud of points shown schematically in 
figure 3. Initially the cloud completely 
fills the region of operation 6. As time 
progresses it changes shape in response to 
the forcing functions and the action of the 
spacecraft control. The maximum projection 
<|>**(t) of this cloud on the cj>-axis at each 
t is the value of the response envelope 
defined by equation (3). In other words, 
the knowledge of the motion of the boundary 
of this cloud is sufficient for the 
computation of the response envelope. 

The motion of this boundary may be 
described mathematically as follows. Let 
the boundary be defined by V(t,x) = 0, with 
V(t,x) < 0 inside the cloud and V(t,x) > 0 outside. Let 

(dV/dt) [t,x,U]_ (t) ,u 2 (t)] denote the time rate of change of V(t,x) along a tra- 
jectory of the system. According to the preceding discussion, the boundary of 
the cloud is characterized by two properties: (i) it is a part of the cloud, 

that is, at each point x on the boundary there is a pair [ui (t) ,u 2 (t) ] such 
that (dV/dt) [t,x,ui (t) ,u 2 (t) ] = 0, and (ii) no trajectory can penetrate it out- 
ward, i.e., at each point x on the boundary and every pair [ui (t) ,u 2 (t) ] in 
Ui(t) x U 2 (t), (dV/dt) [t,x,ui (t) ,u 2 (t) ] < 0. Therefore, V(t,x) satisfies the 
following equation on the boundary 



Figure 3.- Motion of the cloud of states. 


max -nr [t ,x,U! (t) ,u 2 (t) ] = 0 (4) 

Ui(t) x U 2 (t) 

But, the time rate of change of V(t,x) along any trajectory is V t + V x x, and 
x = f [x,Ui (t) ,u 2 (t) ] . Hence, the preceding equation may be expressed as 
follows 


V t + H(t,x,V x ) = 0 (5) 

where 

H(t,x,V x ) = max V x (t ,x)f [x,u x (t) ,u 2 (t) ] (6) 

Ui(t) x U 2 (t) 
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It may be noted that equation (5) is formally a Hamilton-Jacobi equation. The 
boundary of the cloud is the solution of equation (5) such that 
{x: V(0,x) <_ 0} = 0. 

If a closed form solution of the Hamilton-Jacobi equation were generally 
available, the problem of computing the response envelope would be solved, 
because 


cj>**(t) = max <f) (7) 

{x: V (t , x) = 0 } 

Unfortunately, no general solution to Hamilton-Jacobi equation is currently 
available. The best that can be done is to compute points on V(t,x) = 0 by 
means of the canonical equations associated with equation (5) . As already 
noted this results in a lower estimate jf" of the response envelope. Another 
approach is, of course, to try to guess the form of solution of equation (5) 
on the basis of insight in much the same way as one guesses Liapunov functions. 
Since a guess will most likely not be a solution, it is necessary to decide 
what is a meaningful approximation in the present context. From the practical 
point of view, the next best thing to the exact response envelope is a pair 
of lower and upper estimates such that for each t >_ 0, 

<(> _ (t) <_ <|>**(t) <_ cj) + (t). Since may be computed by means of the canonical 

equations, it seems reasonable to consider an approximate solution meaningful 
if, when substituted in equation (7), it yields an upper estimate . A func- 
tion V + (t,x) will be such an approximate solution if the surface V + (t,x) = 0 
encloses the cloud for all t 2 0, that is, if 

{x: V(t,x) <. 0 } C {x: V + (t 3 x)^0} 


for all t >. 0. However, the approximating surface need not be a part of the 
cloud. Hence, V + (t,x) is a solution of the following Liapunov inequality 


vj + H(t,x,vJ <L 0 


(8a) 


with the boundary condition 


(x: V + (0 , x) < 0} 30 


(8b) 


Inequality equation (8a) must hold for all t >_ 0 and all x in 9 such that 
V + (t,x) =0. 

It is very easy to construct functions that solve equation (8) . Simply 
let V + (t,x) = V^( t,x) - V 2 (t) , where V 1 (t,x) is such that for some a < », 

(x: Vi (0,x) <. a} 2 0 (9) 
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and where V 2 (t) is the solution of the following ordinary, first order, 
scalar differential equation with initial condition V 2 (0 ) = a. 


V 2 = max [Vi t + H(t,x,V lx ) ] (10) 

{x e 0: (t , x) = V 2 (t)} 

Then V + (t,x) so defined solves equations (8) as can be tested by direct sub- 
stitution. The corresponding upper estimate may be computed by means of the 
following natural modification of equation (7) . 


cj> + (t) = max cj> (11) 

{x e 0: (t , x) = V 2 (t)} 

Thus, any function Vi(t,x) that can satisfy the boundary condition (9) may 
always be used to compute an upper estimate. Of course, the fidelity with 
which cf> + represents <j>** depends directly on the choice of Vj(t,x). A 
poor choice will result in an overly pessimistic estimate of system 
performance. 


A Prototype of Arbitrary Attitude Control Systems 

In order to construct a Vi(t,x) for an arbitrary three-axis, large-angle 
attitude control system,^, consider the model of such a system in somewhat 
greater detail than shown in figure 1. Associate with every possible motion 
of the system the following time histories: target attitude Ad s , spacecraft 

attitude A^, target angular velocity w^, spacecraft angular velocity w a , 
and spacecraft angular acceleration o_ a . In any practical situation, w^(t), 
w a (t), and o a (t) will be uniformly bounded. That is, if C°(m) is defined as 
the set of piecewise continuous vector functions of time, which are spheri- 
cally bounded by a constant m, there is a triplet (wd ma x> w amax> °amax) of 
constants such that vy e C° (w^max) > ^3. e £ C° ( w amax) > an d 

a a e ^a — ( -'°C cr amax) • This means that any admissible time history Ad s of the 
desired attitude (target) is a solution of the following kinematic equation 
(see eq. (Al) ) . 


A ds = S[w d (t)]A ds (12) 

for some orthogonal A ds (0) and some w^ e £°(wdmax) • Similarly, any possible 
time history of spacecraft attitude is a solution of 


A as = S[w a (t)]A as (13) 

for some orthogonal A as (0) and some e W^. According to equation (A8) , 

equation (2), and the above kinematic equations imply the following expression 
of the time rate of change of system attitude error R. 
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R = S[w a (t) - Rw d (t)]R 


(14) 


This kinematic equation of attitude error will be taken as part of the set E 
of system state equations. Since R is always orthogonal, it could be repre- 
sented using three independent coordinates. This, however, is unnecessary and 
inconvenient at this stage. Instead, all nine elements of R will be con- 
sidered as independent with the stipulation that every R(0) is orthogonal. 

The term -Rwq(t) in equation (14) will be interpreted as the value at time t 
of the forcing function u_i in figure 1. Consequently, U x = C° (wq max ) . 

Because angular acceleration of the spacecraft is bounded, W a has more 
structure than C°(wamax). In fact, any w a is a solution of the dynamic 
equation 


wa = CT a(t) 


(15) 


for some w a (0) <_ w amax and some in the set of admissible angular 

accelerations. Consider any such o^. It is obviously possible to write 
£a = o° + So_, where £° is a specific function and 6o_ = - o°. The inter- 

pretation to be assigned to o° is that it represents that part of the space- 
craft acceleration which helps to reduce attitude error. The difference 6q_ 
will be interpreted as a disturbance. According to the previous discussion 
the magnitude of attitude error R may be represented by the error angle c() , 
and its direction may be represented by the unit eigenvector c of R. Hence, 
any part of angular acceleration that is antiparallel to c may be con- 
sidered helpful. Another helpful part of is one which is antiparallel to 

the angular velocity w a of the spacecraft because it provides damping. For 
these reasons £° will be defined as follows. For every possible motion of 
the system p_° = {(t,z[R(t), w a (t)]): t ^0}, where 


z(R,w a ) = -gi C*Dg2 C ||w a i|)c - g 3 (<l>)w a (16) 

and where gi, g 2 , and g 3 are scalar functions. The perturbation So (t) will 
be represented at every t > 0 by n [t ,R (t) ,w a (t) ]u 2 (t) , where e C°(l), 
and where the function n is chosen so that for every possible motion of the 
system and every t >_ 0 


6 °(t) = - Z [ R 00 >w a (t)] e (n[t,R(t),w a (t)]u 2 (t) : ||u 2 (t)|| <L 1} 


(17) 

The perturbation function n may be a matrix, or simply a scalar. With such 
a representation of spacecraft acceleration £ a , the dynamic equation (15) 
becomes the following: 


w a = z[R(t),w a ] + n[t,R(t),w a ]u 2 (t) 


(18) 
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This equation will be taken as the remaining part of the set E of system 
state equations, and w a will be the remaining part of the state vector x. 
The coordinates of w a are restricted by the condition ||W a || <_ w amax . The 
forcing function u 2 will be identified with the forcing function in the 
disturbance generator as shown in figure 1. Hence U 2 = C°(l). 

The preceding discussion motivates the prototype of attitude control 
systems defined in the following table. 


TABLE I.- PROTOTYPE v/ty ° OF ARBITRARY ATTITUDE CONTROL SYSTEMS 


x = R 12 , X = (R,w a ) 


0 - (x: RR - I , ||w a || <_ w amax ) C 


E =< 


R = S[w a + ui(t)]R 
w = z(x) + n(t,x)u 2 (t) 


Ul = C°(w dmax ), U 2 = C°(l) 


The underlying state space X is 12-dimensional; the state variable is 
denoted mnemonically by x = (R,w a ) . The region of operation 0 is defined 
by the orthogonality condition on the error matrix R and by the spherical 
bound on the angular velocity w a of the spacecraft. Thus, 0 is 
6-dimensional: three for R and three for w a . Clearly, 0 is compact. The 

set E of state equations consists of equations (14) and (18) . The control 
law z(x) is defined by equation (16). It and the perturbation function 
n(t,x) are bounded on 0. The two forcing functions Uq and u 2 belong to the 
sets of piecewise continuous vector functions of time which are spherically 
bounded by w^^ and 1, respectively. The right-hand side of the state 
equations is uniformly bounded on its compact domain. 

In the following very specific sensed 0 is a prototype: it generates 

all time histories of the error angle which are possible for any attitude 
control system having the bounds (wdmax> w amax> a amax) an( ^ the set of 

admissible accelerations for which inequality (17) holds. This inequality can 
be satisfied for any M, however complex, by a proper choice of n(t,x). Con- 
sequently, the response envelope of [wd max > w amax> z ( x ) > n(t,x)] is neces- 
sarily an upper estimate of the response envelope of 

cA= (w dmax , w amax> a amax , Za) for which (17) holds. If n(t,x) is large, the 
estimate will be pessimistic. 
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II 


Explicit Procedure for Computing an Upper Estimate 
for the Prototype 

A very simple V;i(t,x) is now proposed for the computation of an upper 
estimate of the response envelope of sW° by means of equations (9), (10), and 
(11). It is the following. 


Vi(t,x) = [g! (y) + ug 3 (y)sin(t y) ] dy + j J^' Wal1 + y sin(j cf^c 1 ^ 

(19) 

The functions gi, g 2 , and g 3 appear in equation (16). The scalar y is an 
adjustable constant. There is no explicit dependence on time. Hence, = 0. 
The Hamiltonian H(t,x,V lx ) is given by equation (B2) . A Vi(t,x) of the form 
of equation (19) may be thought of as a natural extension to three dimensions 
of the Liapunov function useful for the analysis of one-dimensional servos. 
Thus, the first integral on the right of equation (19) depends only on the mag- 
nitude of attitude error. The second integral depends only on the magnitude 
of angular velocity. The last term represents a coupling which depends upon 
these magnitudes as well as on the angle between the error axis and the 
velocity vector. 

The proposed procedure for computing an upper estimate of the response 
envelope of the prototype -44° is now explicit. It is summarized in the 
following table. 


TABLE II.- COMPUTATION OF AN UPPER ESTIMATE 

Step 1. Obtain Vi using (19) 

Step 2. Choose V 2 (0) to satisfy (9) 
Step 3. Solve (10) for V 2 (t) 

Step 4. Compute $_ + using (11) 


The Hamiltonian H(t,x,Vi x ) needed in step 3 is derived in appendix B. 

Several points regarding this procedure are worth noting. Since Vi 
defined by equation (19) depends only on the triplet [<(), ||w a ||, sin (1/2 c|))c t w a ] 
of scalars, the maximization in step 4 may be performed on a two-dimensional 
surface in a three-dimensional space whose points are related to the states x 
in 0 by the function p(x) = [(f), ||W a ||, sin(l/2 c())c t w a ] . 

The computation in step 3, in general, requires a maximization over a 
five-dimensional surface in 0. If, however, the perturbation function n(t,x) 
depends on x only through the triplet p(x), as above, then this 
maximization can be performed on the same two-dimensional surface as in step 4. 
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Steps 3 and 4 can be performed simultaneously on a digital computer. 
Experience has shown that the computation of an upper estimate in this manner 
takes about one second of computer time. 

Finally, it will be noted that the perturbation enters the computation 
only by way of equation (Bl) needed in step 3. Consequently, when z(x) and 
n(t,x) are selected for any attitude control system^, only that component of 
acceleration need be considered in equation (17) which is in the plane of 

w a and c. 


EXAMPLE 




The following numerical example illustrates the computation and use of r 

upper estimates for the global description of attitude control systems. An 
Orbiting Astrpnomical Observatory (0A0) is considered. Its parameters perti- 
nent to the discussion are summarized in appendix D. The spacecraft is con- 
trolled by means of three identical orthogonally placed motor-reaction-wheel 
combinations . The torque capacity of each motor is &max> and the angular 
momentum storage capacity of each wheel is h ma x. The passive moment of iner- 
tia (moment of inertia of the spacecraft with locked wheels minus the moment 
of inertia of each wheel about its spin axis) is the matrix J a . The maximum 
eigenvalue of J a is jmax- A momentum dumping scheme maintains the total 
angular momentum of the system bounded by h smax . The acceleration is bounded 
by a a max = ^max/jmax, and the velocity is bounded by 

w amax = (hmax - hsmax) /jmax* The torque commanded by the controller on the 
basis of sensor data is J a z(x) , where the control law z (x) is assumed to be: 


z(x) 


1 _ 

2 a amax 


sat(<M s )c 


+ 


1 

w amax 



( 20 ) 


The function sat (4> , <^ s ) = $/$ s for error angle <f> <_ cf> s , and it saturates at 
the value 1 for <j> The saturation angle $ s is a dimensionless combina- 

tion of the dynamic capacities of the spacecraft: 

zn max 

~ 1 i 
Jmax^max 


For the present example of the 0A0, <p s = 0.1 radian. The complete model of 
the system is derived in appendix C. The system has three axes and is 
12-dimensional. Three independent coordinates are needed for the spacecraft 
attitude, three for the desired attitude, three for the spacecraft angular 
velocity, and three for the total angular momentum of the system (spacecraft 
and wheels) . 

The stability of this system has been investigated in reference 4. There 
the system was shown to be asymptotically stable everywhere on its region of 
operation. The plots in figures 2, 3, and 4 of reference 4 show a response of 
the system to a particular initial condition. The questions being considered 


» 


► 
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now are: (1) What is the responsiveness of the system to step changes in the 

desired attitude, and how well does the system follow a desired time varying 
attitude? (2) How significant is the gyroscopic coupling? How sensitive is 
the system to (3) external torque disturbances, (4) variations in system param- 
eters, (5) changes in the form of the control law, and (6) time delays in the 
controller? These questions are resolved by means of upper estimates computed 
on a digital computer using the procedure outlined in table II for various 
appropriate instances of the perturbation function n(t,x). 

Let the system be normalized as follows: time, t t/w amax , angular 

velocity w -* w-w amax . Then, the control law (20) becomes the following. 

z(x) = ' satOMsJc - w a (21) 

This control law is a special case of (16): setting g]_ = sat (<|>,<j> s ) , 

g 2 = 1/4> S , and g 3 = l/<j> s in equation (16) one obtains equation (21), above. 
With these identifications and setting the adjustable parameter y = 1, the 
Vi function defined by equation (19) assumes the following completely explicit 
form. 


Vi(t,x) = p sat (y , < 


>><*♦57 l 1 


COS 


4* s II + sin 


c^w,. 


The results of the computation are summarized below in the same order as the 
questions posed above. 


(1) Responsivenes s of the nominal system . The assumptions for this case 
are: the initial attitude of the spacecraft is arbitrary; the initial angular 

velocity of the spacecraft is arbitrary except that it is spherically bounded 
by w araax ; the initial step change in the desired attitude is arbitrary; the 
desired attitude varies in time arbitrarily except that its angular velocity 
is spherically bounded by wj max ; there are no disturbances, i.e., n(t,x) = 0. 
The result is the family of curves in figure 4. The family parameter b is 
the ratio w dmax/ w amax of the two velocity bounds. The curve b = 0 indi- 
cates the responsiveness of the nominal system to step changes in the desired 

attitude. Thus, for any possible initial condition, and 
regardless of the size and direction of the step change of 
the desired attitude, the attitude error will not be 
greater at any time than indicated by this curve. In fact 
the system is not only globally asymptotically stable, but 
also it is essentially (4> <. 0.01) on target no later than 
5 units of time. Curves with b > 0 indicate how well the 
nominal system follows the desired attitude when that is 
varying in time. Thus, the curve b = 0.2, for example, 

, t eS , shows that for any possible initial condition, and regard- 

0 2 4 tw 5^ 8 10 less of how the desired attitude varies in time, so long 

as its angular velocity remains spherically bounded by 
Figure 4.- Nominal system. 0.2 w a jjj ax , the attitude error will not be greater than 

that indicated by the curve. It is emphasized that a 
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curve in figure 4 is not a response of the system. Rather, it is a descrip- 
tion of all (there are infinitely many) responses of the multivariable, 
multidimensional, nonlinear, and nonautonomous system with infinite sets of 
initial conditions and forcing functions. 

(2) Significanc e of the gyroscopic coupling. It is assumed that initial 
spacecraft attitude is arbitrary and its angular velocity is spherically 
bounded by w amax . The desired attitude is stationary after an arbitrary step 
change. The system is preloaded with a total angular momentum which is spheri- 
cally bounded by h smax but is otherwise arbitrary. There are no net exter- 
nal torques. The appropriate perturbation term is given by equation (C6) . 

The result is the family of curves in figure 5. The family parameter b is 

the ratio h sniax /h max of the bound on the total angular 
momentum to the momentum capacity of the reaction wheels. 
The curve b = 0 corresponds to the case of no gyroscopic 
coupling. The curve b = 0.3 indicates that for the case 
of the 0A0 considered, gyroscopic coupling is not very 
significant even when the system is loaded with as much 
as 30 percent of its angular momentum storage capacity. 

(3) The effec ts of external torque disturbances. 

The assumptions in this case are as in case (2) except 
0 2 4 6 8 1 o that net external torque is allowed to be nonzero. The 

„ . only restrictions are that this torque is spherically 

effects. bounded by &smax> and that it does not overpower the 

momentum dumping system which maintains the total angular 
momentum bounded by h smax . The appropriate perturbation 
term is given by equation (C 7 ) . The result is the pair 
of curves in figure 6. The lower curve is for 
hsmax = ^smax = 0. The upper curve is for 
hsmax = 0.3 h max and f smax = 0.1 £ m ax- Thus, for any pos- 
sible initial condition, attitude error will be less than 
about 10° after 6 units of time despite external torques 
which may be as large as 10 percent of the torque 
capacity of the controller. 

02468 |0 - , _ „ ... .. 

tw mox — (4) Sensiti vity to variations m system parameters. 

Figure 6.- External torque Until now the spacecraft was assumed to be controlled 

disturbances. exactly by the control law (20). That is, it was assumed 

that the motor-reaction-wheel combinations were placed exactly on the orthogo- 
nal body axes, that the error angle <f>, the error axis c, and the angular 
velocity w a of the spacecraft were computed from sensor data without error, 
and that the matrix J a was known exactly. Obviously, in practice, error 
will be present because motor-reaction-wheels will be slightly off axis, sen- 
sor characteristics will vary, signal and power amplifiers will drift, the com- 
puter will have roundoff errors, and the moment of inertia of the spacecraft 
will vary because of shifting mass caused by thermal distortion and movement 
onboard of instruments and personnel. The purpose of the present subsection 
is to investigate the effects of such errors. 

Let the perturbation in the active parameters be denoted by the column 
6p of appropriate dimension, and let the set of all possible 6 p be P 
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which may be assumed to be infinite. The actual torque generated by the con- 
troller is a function, say J°z(x,6p), of the state x and perturbation 6p. 
Then the effects of perturbation on the angular velocity of the spacecraft may 
be described by the function m(x,6p) = j£ 1 J^[z(x,0) - z(x,6p)], where J a is 
exact moment of inertia of the system. For a fixed state x, let the set 
M(x) = [m(x,6p) such that 6p e P] and let the set N(x) = [n(t,x)u 2 (t) such 
that ||u 2 (t)|| <. 1] where n(t,x) is chosen so that for t > 0 M(x)C N(x) for 
all x in 0, the region of operation of the system. Clearly, regardless of 
how 5p varies in time, its effect on the angular acceleration may be 
accounted for by the perturbation term n(t,x)u 2 (t), where u 2 e C° (1). The 
significance of the perturbation in system parameters may now be estimated 
with an upper estimate of the response envelope. Two examples of 
perturbations are given next. 



Figure 7 • “ Spherical 
error in acceleration. 


(i) Spherical perturbation in angular acceleration. 
n(t,x) = biiz(x,0)ll. The constant b is the intensity of 
perturbation, and z(x,0) is the nominal control law z(x) 
given by equation (20) . The curves in figure 7 are upper 
estimates for various intensities b. The plots indicate 
that spherical errors in acceleration of the order of 
10 percent affect the performance of the system little. 
This means, for example, that 10 percent drift in moment 
of inertia, motor and power amplifier gains, or a mis- 
alinement of the motor-reaction-wheels of about 3° is not 
very detrimental to the performance of the system. Even 
when such errors are large enough to cause 30 percent 
errors in acceleration, the system remains globally asymp- 
totically stable, and it is essentially on target after 
about 8 units of time for any possible initial condition. 



(ii) Spherical error in the error axis. Suppose that 
the difference between the true error axis c and its com- 
puted approximation c° is spherically bounded by b, 
that is, suppose that ||c - c°|| <_b. Then the correspond- 
ing n = bsat (cj) , cj) s ) /cj) s . Figure 8 shows the computed upper 
estimates for several intensities b. The plots indicate 
that the system is not very sensitive to errors in the 
error axis. In fact, even when these errors are as large 
as 40 percent, the system remains asymptotically stable 
everywhere on its region of operation, and it is essen- 
tially on target for any possible initial condition after 
about 8 units of time. 


Figure 8.- Spherical 

error in error axis. (5) Sensitivity to changes in form of control law . 

In the preceding subsection perturbations in angular acceleration were assumed 
unavoidable and largely undesirable. Another point of view is possible. Sup- 
pose that all system parameters are essentially fixed, but that it is desir- 
able for reasons of simplicity of implementation to use a control law which is 
different from that given by equation (20) . In such a case an upper estimate 
to the response envelope of the modified system may be computed as in the 
preceding subsection if the given control law is considered a sum of 
equation (20) and a perturbation. The next two examples illustrate this point. 
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(i) A processor of signals from attitude sensors. Let the attitude of 
the spacecraft be measured with some arbitrary combination of star trackers, 
inertial gyros, sun sensors, etc. Let the outputs of all sensors be arranged 
in a single column matrix y of an appropriate dimension, and let the com- 
plete package of sensors be described by y = g(A as ,y°) where A as is the 
attitude of the spacecraft and y° is the set of inertial coordinates of the 
sensor targets (i.e., guide stars, inertial directions, sun, etc.). Let y, 
A ds , and y° be combined in a processor whose output is some given function 
f (y,Ads >y°) • Finally, suppose that this output is used instead of sat (<j> ,<j> s )c 
in equation (21); that is, suppose that the new control law is the following 
modification of equation (21) . 


Z * (X) = “ *7 [ £(y ’ A ds>y°) + w a ] 

Then the perturbation in z (x) is z (x) - z*(x) = [f - sat (4>,4> s )c]/4> s . Con- 
sequently, the upper estimates given in figure 8 describe this modified system 
for constant A ds and y°, and, of course, <p s = 0.1. The intensity 

£ [gCRA ds ,y°),A ds ,y°] - sat (4>,<P s )c 
“sat (^s) 

Thus, one may conclude, for example, that if reaching the target in 8 units of 
time is sufficiently fast for the mission of the spacecraft, any combination 
of attitude sensors and processor with b <_ 0.4 may be considered acceptable. 
Two possible reasons one might be interested in a case with b > 0 are (1) sim- 
plification of the control law, and (2) determination of effects of failure of 
part of the output of the sensor package. 

(ii) Control with control moment gyros. Suppose that the spacecraft is 
to be controlled not with reaction wheels but, rather, with a set of control 
moment gyros. Let the active gimbal angles of all the gyros in the package be 
arranged in a column y of an appropriate dimension, and let the spacecraft 
coordinates of the total angular momentum of all gyros be denoted by h^. 
Assuming that the total angular momentum of each gyro may be adequately approx- 
imated by its spin momentum, we may express h a as a function of y. Let 
h a = h(y), where h(y) is one-to-one on Y, and h(Y) is a solid sphere with 
radius h max . Then the control torque is -hyy. Let the gyro gimbals be 
driven through a processor so that y = -F(y)J a z(x); where F (y) describes the 
processor, J a is the moment of inertia of the spacecraft, and z (x) is the 
control law given by equation (21) . Then the new control law is given by 
z*(x,y) = J a x hyF (y) J a z (x) , and the perturbation 

z(x) - z*(x,y) = [I - J a x hyF (y) J a ] z (x) . Consequently, the upper estimates 
given in figure 7 describe the spacecraft controlled with a set of control 
moment gyros when <j> s = 0.1. The intensity is 


b = 


max 

r» t 


b = max III - Ja 1 h y F(y)J a ll 
Y 
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Two possible reasons one may be interested in a case with b > 0 are 

(1) simplification of the control law, and (2) determination of effects of 

failure of a number of gyros in the package. 

(6) Sensitivity to t ime delays . It is assumed that the angular velocity 
of the spacecraft is not measured directly but, rather, that it is computed 
from attitude data supplied by attitude sensors, and that this computation 
yields the exact angular velocity delayed by A units of time. The resulting 
control law is the following modification of equation (21) . 

z* (x. A) = - [sat(4>,<)> s )c + w a (t - A)] 

The perturbation z(x) - z*(x,A) = [w a (t - A) - w a (t)]/t)> s . But, 

w a (t - A) - w a (t) = -J t t+A z*[x(t),A]dt 


Consequently, the effects of small time delays are described by the curves in 
figure 7, where the intensity b = A/<j) s = 10 A. In particular, the curve 
b = 0.1 describes the 0A0 with a real time delay of 3.33 seconds. An estimate 
not restricted to small A may be computed using 


n = f-^-lmax z*(x,A) 

+s e 

CONCLUSION 


The purpose of an attitude control system is to force the spacecraft to 
track a target attitude regardless of disturbances and initial conditions. 
Therefore, there is always the problem of describing how responsive the pro- 
posed system is to inputs, how unresponsive it is to disturbances, and how 
quickly it overcomes the initial conditions. The present note proposes the 
concept of response envelope as a solution to this problem. In addition, 
methods for computing lower and upper estimates of the response envelope are 
presented. These methods are useful because, although the response envelope 
is theoretically simple, it is difficult to compute in practice. The lower 
estimate is computed by means of the standard theory of functional maximiza- 
tion. The upper estimate is computed by means of a Liapunov inequality and a 
prototype of attitude control systems. A proposed attitude control system may 
be judged unacceptable on the basis of a lower estimate. Conversely, it may 
be judged acceptable on the basis of an upper estimate. The numerical example 
presented in the note shows that the response envelope and its estimates are 
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useful for describing the responsiveness of complex systems whose behavior 
cannot be ascertained by a direct enumeration of cases. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif. 94035, June 20, 1968 
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APPENDIX A 


SOME USEFUL PROPERTIES OF THREE-DIMENSIONAL ROTATIONS 


(1) Let A be a 3 x 3 orthogonal matrix whose determinant is +1. It can 
be interpreted as a rotation of an orthonormal triplet a relative to an 
orthonormal triplet s. Then, if x s and x a are the coordinates of an 
arbitrary vector x relative to the s and a triplet, respectively, x a = Ax s . 

(2) Suppose that A is a function A(t) of a real variable t. Then it 
can be shown (see ref. 4) that there is a column matrix u(x) such that 


S(x) = 


dA _ 
dx 

S[u(x)]A 


(Al) 

t X, 

the skew 

symmetric matrix 

0 

*3 

-X 2 \ 


-*3 

i 

0 

X 1 

(A2) 

\ x 2 

" x l 

° y 



If x is interpreted as time then u(x) gives the a-coordinates of the 
angular velocity of a relative to s. 

(3) According to Euler ! s theorem on rotations, A may be interpreted as 
the result of a rotation from identity about a fixed axis. Hence, A is the 
solution of 


dA 

dx 


S (c) A , 


A (0) = I 


at some x = cf>, and some constant c such that c^c = 1; <j> will be referred 
to as the angle of A, and the eigenvector c as its axis. The solution is 

A = = I + sin <j)S(c) + (1 - cos 4>)S 2 (c) (A3) 


The angle (j) and the axis c may be computed from the elements of A by the 
following formulas: 

cj)(A) = arc cos [tr(A) - 1]\ (A4) 

[0,tt] 12 J 
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c(A) = y cosec <j> 


_a 32 \ 


I a 2 3 
a 31 _a l 3 I 

a 12 _a 2 l / 


(A5) 


It may be noted that ^(At) = <j>(A), cfA 1 ") = -c(A), and, of course, Ac = c. 

(4) Let 0^ be the angle between the ith vector of a and the ith 
vector of s. Then 0i <_ 4 > (A) . Indeed, 


cos 0i = cos <j>(A) + [1 - cos <t>(A)]c|CA) 


Hence, cos 0i >. cos <t>(A). 

(5) Consider all paths from I to A. Each satisfies the differential 
equation 


dA 

dx 


S[u(x)]A 


with A(0) = I, A (x i ) = A, and some function u. Then 

<j) (A) <_ /J 1 ||u(T)lldx 


so that <j) (A) may be considered to be the minimum angular distance between the 
a and the s triplet. This may be shown as follows. The Hamiltonian for 
the minimization problem is 


H = tr{PtS[u(T)]A} + p 0 ||u(x)|| 


and 


P 0 = 0 

P = S [u (x) ] P 


Thus, A and P have the same transition matrix $(t) . Hence, 


H = 2u -t (x)$ (x)k + p 0 ||u(x)|| 
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where k is constant. The direction of u(x) which minimizes H is 


u(t) _ -2$ (x)k 
I|u(t)|| " p 0 

But this means that the direction of u(x) is fixed in the s-triplet. 
Therefore, u(x) is at each x the eigenvector of A(x), and the conclusion 
follows . 

(6) For any three-dimensional rotations A and B, 


<j> (AB t ) <_ <j>(A) + 4> (B) 


(A6) 


Suppose the contrary, and let ABt = c and B = D^. Then it would be true that 
4> (C) > 4>(D) + That is, the angle of the composite rotation: from I 

to D, followed by D to C is smaller than the angle of direct rotation 
from I to C. This, according to the preceding section, is impossible. 

(7) Consider the set of all three-dimensional rotations. For any A and 
B in this set let the following function be defined 


<f> (B,A) = <f>(ABt) 


(A7) 


The function c()(B,A) so defined is a metric on the space of three-dimensional 
rotations. Indeed, (i) cf)(B,A) is positive; (ii) (})(B,A) = 0 if and only if 
A = B; (iii) <J)(A,B) = <f>(B,A); (iv) <J> CB, A) + <|>(A,C) <J> (B, C) . The triangle 

inequality follows from the preceding section: 


<j>(B,C) = <f> (CB 1 ) = <j> (CA t AB t ) = [CA 1 " (BA 1 ) <_ ^(CA 11 ) + 4> (BA 1 ) = <1>(A,C) + <f>(B,A) 


(8) Let A = A(t) and B = B(x), where x is real. According to prop- 
erty (1) u and v are such that A = S(u)A and B = S(v)B. Let C = AB^. Then 

C =ABt + AB^ = s (u) ABt -ABtS(v) =S(u)C -CS(v) =S(u)C -CS(v)CtC = S(u)C- S(Cv)C 


Hence 


dC 

dx 


= S (u - Cv)C 


(A8) 


(9) Let A = e4 ) S(c) j and dA/dx = S(u)A. Then it follows (e.g., 
appendix B of ref. 4) that 
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is the vector part of the quaternion of A. 


CA9) 

(A10) 
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APPENDIX B 


THE FORM OF THE HAMILTONIAN 


The projection of state point velocity on the gradient Vi x is Vi x x. 
The Hamiltonian appearing in equation (10) is according to equation (6) the 
maximum of V lx x over Ui (t) x U 2 (t). For defined by equation (19), 


Vix* = [giO) + PgsWsinft <(>)] 


dc|) 

dt 


2g 2 (Hw a ll) 


_4_ 

dt 


nw a ir + p 


_d_ 

dt 


sin ( 


tf.jcV 


where the derivatives are evaluated along the trajectory. They may be 
obtained as follows. 

(i) According to equation (A9) and the kinematic equation in table I, 


dcj> 

dt 


c t [w a + u x Ct)] 


(ii) Since 


Iw a ll 2 = W^w a , 



2w^w o 
a a 


(iii) From equation (A10) , it follows that 
dT [ sin ^ ^ ctw a] = w a A [ Sin ^ +) c ] + Sin ^ ^"a 

= w a j (sin(~ <j>) S [w a + Ul (t)]c + cos (j c|>) [w a + u x (t)] j 

+ sin <p)c\ 


But w a = o & (t). Hence, 



Vix* = [y sin(j*)c + g 2 + [«i(+) + Mg 3 (+)sin(j +) ] c l w a 

+ j y cos(j (j)) II w a tl 2 + ui(t) [giC+Dc + pg 3 (<f>)sin(j cf>)c 
+ j y sin (y 4>)sCc)w a j 

If c? a (t) = z[R(t),w a (t)] + n[t,R(t) ,w a (t)]u 2 (t) , and z(R,w a ) is given by 
equation (16) , 

Vixx = -y gi WgaCll w a ll) sin (^ *) -[ g2 o?wJ) ~ I y cos (j*)] l|Wal ' 2 

+ ui(t)[gi(*)c + y g 3 (<f))sin(y *Jc + \ y cos (j <fr)w a + j y sin(j *JS(c)w a J 

+ (t^n 1 (t,R,w a ) [ + sin(t *)c] (Bl) 

The Hamiltonian in equation (10) is 


H(tjX,V]!x) - max ^i x x 

C ( w dmax) x C°(l) 


(B2) 


where Vi x x is given by (Bl), above. 

For the particular control law (21) used in the example. 


V ix x = - j- sat(<£,<l> s )sin(y <p] - [ 1 - j cos (j <f>) 


II w„ 


u\ (t) f sat (4> , 4> s ) c + -t sin(y <j>)c + j cos (t <t>) w & + j sin(j <|>) S (c)w & J 


+ U2(t)n t (t,R,w a ) 


f*s w a + sin 


in(| <J>)c j 
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APPENDIX C 


DETAILED MODEL OF SYSTEM USED IN THE EXAMPLE 


The dynamic equation of a spacecraft controlled by means of an arbitrary 
angular momentum exchange and storage device (i.e., reaction wheels, control 
moment gyros) may be obtained as follows. The total angular momentum of the 
system is the sum of the angular momentum stored in the main body and that 
stored in the controlling device. Denote the inertial coordinates of the 
total angular momentum by h s , the spacecraft coordinates of the portion 
stored in the device by h^, and the spacecraft coordinates of the moment of 
inertia of the main body by J a . Then, 


^as^s ^a w a + ^a 


(Cl) 


Taking the time derivative, using equation (13), and rearranging terms, one 
obtains the following equation. 



+ J a * S (w a ) A as h s + J a i A as h s - J a j J a w a 


- 1 / 


r-l 


(C2) 


Thus, the angular acceleration of the spacecraft is the sum of four terms. 
The first will be interpreted as the control acceleration, and the control 
torque will be defined as follows. 



(C3) 


The second term will be interpreted as gyroscopic coupling. The third is 
caused by external torque. The body coordinates of the external torque will 
be denoted by £|; hence. 


^a A as^s 


CC4) 


The fourth term is present only when the moment of inertia of the main body 
varies in time. 

Any practical angular momentum exchange and storage device is limited in 
both exchange rate and storage capacity. Thus, the torque of motors driving 
reaction wheels is limited as is the speed of the wheels. Similarly, the 
rates with which gyro gimbals can be driven are limited by available torque 
capacity, and the geometry defining the arrangement of gyros in the package 
imposes a limit on angular momentum storage. This fact will be accounted for 
by the following limits: || £ a || <. & max and l|h a l| <_ h max . In addition, it will 
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be assumed that angular momentum is dumped in such a way that II fcfll <_ 5-emax 
llh s ll <. h S max- Equation (Cl) implies that II h a ll <_ h max if 

II w a || <. wamax = (h max - h smax )/j ma x, where j max is the maximum eigenvalue of 

J a . 

Equations (2), (12), (14), (C2) , (C3) , and (C4) may now be combined as 
follows : 

TABLE III.- DETAILED MODEL OF THE SYSTEM 
X = ft 24 , x = (A ds ,R,w a ,h s ) 

0 = (x: A ds A ds = I, RR = I, ||w a || <_ w amax , II h s ll <. h smax ) C & 12 


E = 


Ads - S [wd (t) ] Ads 


R = S[w a - Rw d (t)]R 


w„ 


= J"l(t)£?(t) + J a X (t) 


S(w a )RA ds h s + A®(t) + J a (t)w a 




Hd ^CC 0 ^), £ e L c CC»(t lax ), £ c I e CC»« eiaI ), J, S JCC‘, 
and L c x L e x J is such that any x(t) e 0 for all t >_ 0 . 


The underlying state space of this model is 24-dimensional. However, the 
region of operation 0 is 12-dimensional. Three dimensions are for target 
attitude A ds , three for attitude error R, three for angular velocity of the 
spacecraft w a , and three are for the total angular momentum of the system h s . 
The motion of the system is given exactly by this model for any admissible ini- 
tial condition and time variations of target velocity, control torque, exter- 
nal torque, and moment of inertia. It may be noted that in the absence of 
external torque, the total angular momentum is a constant of the motion, and 
the last state equation may be dropped. If, in addition, the total angular 
momentum is zero, then the gyroscopic term is absent. 

In the example given in the main text, J a 1 (t))l a (t) is assumed to be of 
the form z(R,w a ) given by equation (16). Hence, the perturbation appearing 
on the right-hand side of the inequality (17) is the following. 


6a(t) = J- : (t) [s(w a )RA ds h s + £|(t) + J a (t)w a 


(CS) 
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Consider just the gyroscopic term. It is, clearly, bounded by 
^smax^max^' w all • This bound could be used as the perturbation function n 
in figure 4, The perturbation 60 (t) = n(t,x)u 2 (t) enters in the computation 
of an upper estimate only by way of equation (Bl) . Therefore only the compo- 
nent of 6a(t), which is in the plane of w a and c, is significant. For this 
reason if the scalar n is replaced by the matrix N, defined as follows, a 
finer estimate will be obtained. 


"tv ■ i [ s( " a) * a "” a " 1 ] CC6) 

Jmax^ 1 + a J L J 

where a = (jmax/jmin) - 1- This may be justified as follows. For any 
(A ds , R , w a, h s) in e > 

J a ^ S (w a ) R Ad s h s e ( y : y - hsmax^a^Cwg^ A anc ^ -^-1} 


But that set is included in 


y: y = 


JI smax 
3 max 


S(w a )Ai + a|| w a || A ^ 


and || A ill <_ 1 and l|A 2 ll <. 1 


This set, in turn, is included in 


{y: y = N(w a )A and ||A|| <_ 1} 


Thus, the perturbation N(w a )u 2 (t) generates all possible cases of the 
gyroscopic term. 

The effects of external torque as well as the gyroscopic coupling can be 
represented by the following perturbation. 

6a (t) = N(w a )u 2 (t) + u 3 (t) (C7) 

Jmm 

where both u 2 and U3 belong to C°(l). 


29 


APPENDIX D 


PARAMETERS OF THE SYSTEM USED IN THE EXAMPLE 


maximum principal moment of inertia j max 
minimum principal moment of inertia jjnin 
angular momentum storage capacity ^max 
maximum torque capacity ^max 


4>s 


2h. 


max 


Jmax^max 


0.1 


1.4xi0 3 kg-m 2 
l.OxlO 3 kg-m 2 
4.2 N-m-sec 
0.25 N-m 


a 


3 max ~ 3 min 
-’min 


0.4 


w 


o 

max 


hmax 
3 max 


3 milliradians/sec 
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